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Abstract: An improved method of energy calibration of position-sensitive silicon detector is presented. Instead of 
a parabolic function used in the traditional method, a new function describing the relation of position and energy is 
introduced and better energy resolution is achieved. Using the new method, the energy resolution of the 8.088 MeV a 
decay of ^^^Rn is determined to be about 87 keV (FWHM), which is better than the result of the traditional method, 
104 keV (FWHM). In addition, different functions can be tried in the new method, which makes the calibration of 
detectors with various performances possible. 

Key words: position-sensitive silicon detector, energy calibration, correction of linear calibration, energy resolution 
PACS: 29.40.Gx, 29.30.Ep, 25.70.Hi 


1 Introduction 

Fusion-evaporation reactions are generally used in the 
synthesis of heavy nuclei. Evaporation residues are sepa¬ 
rated from primary beam and products of transfer reac¬ 
tion by the electromagnetic separation device [T] . Subse¬ 
quently they are implanted into position-sensitive silicon 
detector. New isotopes are usually identified by estab¬ 
lishing a decay chains leading to the known transitions. 
The position-sensitive silicon detectors (PSSD) are ex¬ 
tensively applied in the experimental studies of heavy nu¬ 
clei, where position and time correlation measurements 
between implanted nuclei and subsequent a decays are 
performed. The PSSD used in the present work is a type 
of Xl-300 manufactured by Micron corporation (UK). 
This detector is ion-implanted with an active area of 50 
mm X 50 mm and a thickness of 300 fim. The silicon 
detector is fully depleted when a reverse bias of about 30 
V is applied. One surface of the detector, with a well- 
distributed resistive layer on it, is divided into 16 strips, 
3.125 mm wide for each strip. Signals are output from 
the two ends of each strip. The other side of the detec¬ 
tor evaporated with a thin layer of Aluminum is used to 
get the total signal of the silicon. In the present work, 


three PSSD detectors are installed side by side to form 
an implantation detector. 

In the following two sections, a new method of en¬ 
ergy calibration of PSSD will be introduced, and then 
the calibration results of the new method and traditional 
method mm will be compared through some long-lived 
a emitters produced in the ^°Ne+^°®Bi reaction. 

2 Energy and position calibration of 
PSSD 

Signals output from the two ends of each strip are 
amplified in preamplifier and main amplifier before in¬ 
put into ADC (Analog to Digital Converter). The po¬ 
sition of an event along the strip is determined by the 
resistive charge division [4], while the resolution of posi¬ 
tion depends on the performance of detector and noise 
from electronics. The energy registered in one strip is 
the sum of energy output from the two ends of the strip. 
The position perpendicular to the strip is given by the 
strip number. 

In Fig. I, two-dimensional spectrum of signals from 
the two ends of a strip for external a sources is plotted. 
The sources are: (5.157 MeV), ^^^Am (5.486 MeV) 
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and (5.805 MeV) [5], which correspond to three 

groups of events in the middle of Fig. 1. 

Three abnormal phenomena are observed in Fig. 1. 
At first, each of the three groups of events is not lin¬ 
ear distributed but curve. This is mainly caused by the 
mismatch between input impedance of preamplifier and 
resistance of our detector. Events in the center of strip 
have much longer rise time than those in the two ends of 
strip IMS]. As a result, the ballistic deficit is much more 
serious for the events in the center of strip. Addition¬ 
ally, some other reasons may also lead to this distortion 
[9] . This distortion can be reduced by an appropriate de¬ 
sign of preamplifier. Secondly, the deviation caused by 
noise is position dependent. Three types of noise are pre¬ 
sented for a similar detector in reference [ 10 ], which are 
the noise from the input transistor of amplifier, thermal 
noise from the resistive layer and shot noise produced 
by the leakage current of the detector which is mainly 
caused by radiation damage. The amplifier noise and 
shot noise are larger at the ends of the strip than at the 
center of it nni. At last, the data points are limited to a 
finite sector region in Fig. 1, but not extend to the co¬ 
ordinate axes. It results from the finite input impedance 
of preamplifier and can be restored by a simple relation 
m- This situation can be improved by an appropriate 
design of preamplifier and the increase of integral time of 
main amplifier. However, the increasing of integral time 
may result in worsening accuracy of position measure¬ 
ment. A simple model is developed to simulate these 
phenomena. The simulated results can reproduce the 
experimental data curve. 

Due to the first distortion mentioned above, a linear 
calibration of energy is not enough. In order to improve 
the result of linear calibration, parabolic fitting is pro¬ 
posed in the traditional method, in which the systematic 
variations of vertex and curvature of the parabola with 
a energy must be considered W- A new method of the 
improvement, which could be applied more universally to 
silicon detectors based on the resistive charge division, is 
presented below. 


2.1 Linear calibration of energy 


The energies Ei and E 2 extracted from the two ends 
of a strip are proportional to their ADC amplitudes Xi 
and X 2 , respectively. 


El =aiXi + 6i, 


E 2 — <^2^2 T ^2- 

The total energy is 


E — El + E 2 . 


Then, a relation of E and Xi, X 2 can be obtained 


X2 = - Xi-\ - E- 

a2 a2 


bi + 62 
a2 


(1) 

( 2 ) 

( 3 ) 

( 4 ) 



Fig. I. A representative two dimensional spec¬ 
trum of the raw energy signals (ADC amplitudes) 
from the two ends of one strip. 


Three parameters /co=-ai/a 2 , 1 / 0.2 and bsum={bi-\-b 2)/02 
can be determined by linear fitting of data points (Fig. 
I), where bi and 62 are parameters to be determined 
later. Then, the energy obtained by linear calibration 
can be written as 


E = a2{-koXi + X 2 + bsum)‘ (5) 

It is supposed that the center of noise distribution 
(shown in the lower corner of Fig. I), namely the maxi¬ 
mum of noise distribution density, is situated at the mid¬ 
dle of a strip [2| . This assumption is confirmed by a test 
in reference m- So the output energies Ei and E 2 cor¬ 
responding to this center are equal. From Eqs (I) and 
( 2 ), we have 


aiX^i+61 =a2X^2 + ^2, (6) 

where x^i and Xn 2 are the ADC amplitudes of noise at 
this center, then bi and 62 can be obtained from bsum 
and Eq. ( 6 ). 

2.2 Position calibration 

As in reference [3], the position along a strip can be 
determined by the difference between energies Ei and E 2 
extracted from the two ends of the strip. The relative 
position T] is given as 


El — E2 
El + E2 


( 7 ) 


where Ei and E 2 are the results of linear calibration of 
energy, shown in Eqs. (I) and (2). Relative position 
gives the position relative to the geometrical center of a 
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strip, the range of which is between -1 and 1. Introduce 
Eqs. (1) and ( 2 ) into Eq. (7), we have 

_ aiXi-a2X2 + {bi-b2) . . 

CLiXi-\- CI2X2 “h (^1 “h ^2) 

The relation between real position and relative position 
is 

P=^v, (9) 

where L is the length of strip, L = 50 mm here. 

2.3 The improvement of linear calibration by 
new method 

The energy resolution can be improved by a new 
method based on the result of linear calibration. 

The abnormal phenomena presented above can still 
be observed in Eig. 2, where the relation of relative po¬ 
sition 77 and energy E obtained from linear calibration 
is described. The events close to the ends of each strip 
(see Eig. 2) cant be distinguished for different a sources, 
therefore only the events around the central part, encir¬ 
cled with colored lines (see Eig. 2), are chosen in the 
improvement of liner calibration. When an event occurs 
at one end of a strip, the output at the other end will 
be zero ideally, the range of 77 will be 1 or - 1 . How¬ 
ever, as mentioned above, due to the existence of input 
impedance of preamplifier, the range of relative position 
obtained by the resistive charge division is narrowed. It 
can be seen in Eig. 2 that the relative position 77 spans 


approximately between -0.5 and 0.5. 
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Fig. 2. Two dimensional spectrum of relative 
position 77 and E, the energy obtained from lin¬ 
ear calibration. Three groups of events circled 
by lines of different colors define good events of 
the three external a sources, which can be distin¬ 
guished from each other. 



Fig. 3. Two dimensional spectrum of g and E/El 
for the selected events in Fig. 2, see text for ex¬ 
planation. 


From Eqs. (3) and (7), we have 



E 2 

Ex 1 + ?? ’ 

( 10 ) 

make 

then 

Z 

5=1, , 

l + r? 

E 

( 11 ) 


'Ex^^' 

( 12 ) 


Eq. (12) indicates a linear relation between EjEi and g. 
But the data deviate from this simple relation. 



Fig. 4. Two dimensional spectrum of relative po¬ 
sition 77 and E, the energy improved by the new 
method. 
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In Fig. 3, the selected events of the three a sources 
in Fig. 2 are plotted in the g — EjEi coordinate system. 
These three groups of events nearly overlap and follow 
a smooth curve. As a result, a new function, f{g), is 
introduced 

^ = /(5) (13) 

to describe the actual relation between E/Ei and g. 

Different functions have been tried fitting data points 
in Fig. 3, and a function with the following simple form 

f{g) = Ag^ + Bg + Cg^'^ + D (14) 

is found to be among the best. The results of linear cal¬ 
ibration are improved by this function and new results 
are shown in Fig. 4 for all the data in Fig. 2. 

There is some energy loss in the dead layer for the ex¬ 
ternal a sources due to the existence of dead layer of the 
detector. So, for the internal a sources, the energy cali¬ 
brated by the external a sources is larger than their real 
value. In order to reduce this difference, some long-lived 
residues in the ^°Ne+^°^Bi reaction are used to make a 
linear correction for the result obtained above. We have 

E = [Ei{Ag^ + Bg + Cg^^^ + D)]K + M (15) 

where K and M are the parameters of this linear correc¬ 
tion. 

3 Comparison of the calibration results 
obtained by new method and tradi¬ 
tional method 

This new method has been applied in the experimen¬ 
tal study of the ^°Ne+^°^Bi reaction performed on the 
gas filled separator SHANS in Lanzhou [12]. The beam 
energies were 111.4 MeV and 122.6 MeV for ^°Ne. 



Fig. 5. A comparison between calibrated results 
by new method (red line) and traditional method 
(dashed blue line). 


Table 1. Comparison of energy resolutions 


nuclei 

a decay energy/MeV 


FWHM/keV 

new traditional 

2iiAt 

5.870 

77 

79 

212Rn 

6.264 

78 

80 

213pj. 

6.775 

80 

85 

211po 

7.450 

83 

94 

212At 

7.679 

88 

98 

212 m At 

7.837 

99 

114 


8.088 

87 

104 


The calibrated results by the new method and tradi¬ 
tional method for the a decays of some isotopes produced 
in the ^°Ne+^°^Bi reaction are presented in Fig. 5, while 
the energy resolutions obtained are listed in Table 1. It 
is obvious that the new method produces better results 
than the traditional one, especially for the a decays with 
higher energies. 

For the a decay (8.088 MeV) of the transfer reac¬ 
tion residue ^^^Rn in the ^°Ne+^°®Bi reaction, the en¬ 
ergy resolution is determined to be 87 keV (FWHM) 
which is better than the result of traditional method, 
104 keV. The position resolution is determined to be less 
than 2 mm (FWHM), similar to the result of traditional 
method. The position resolution is determined from the 
differences between the positions of short-lived implanted 
nuclei and those of the subsequent a decays. The res¬ 
olution of position is not improved comparing with the 
result of traditional method, which may result from the 
big noise(seen in Fig. 1). Noise-reduction can improve 
the resolutions of energy and position. 

4 Summary 

A new method of energy calibration of position- 
sensitive silicon detector is introduced based on the ob¬ 
servation that the distributions of data in the g — E/Ei 
coordinate system are almost identical for different a 
sources (see Fig. 3). Compared with traditional method, 
the procedure is simpler and the obtained energy reso¬ 
lutions are better. Different functions can be tried in 
this new method making the calibration, therefore de¬ 
tectors with different performances can be calibrated by 
this new method, especially for detectors, the character¬ 
istic (namely the two dimensional spectrum of position 
and energy) of which violates parabolic shapes. 
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